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ABSTRACT 

If the first stars formed soon after decoupling of baryons from the thermal cosmic background 
radiation (the CBR) the radiation may have been last scattered in a cloudy plasma. We discuss 
the resulting small-scale anisotropy of the CBR in the limit where the plasma clouds are small 
compared to the mean distance between clouds along a line of sight. This complements the 
perturbative analysis valid for mildly nonlinear departures from homogeneity at last scattering. 
We conclude that reasonable choices for the cloud parameters imply CBR anisotropy consistent 
with the present experimental limits, in agreement with the perturbative approach. This means 
the remarkable isotropy of the CBR need not contradict the early small-scale structure formation 
predicted in some cosmogonies. 



Subject headings: cosmology: cosmic microwave background — galaxies: formation 



1. Introduction 

Cosmogonies in which the baryons were concentrated in clouds at the epoch of last scattering of the 



thermal cosmic background radiation (the CBR), as in isocurvature models ( Peebles 1994 , Peebles 1997a ) 
are in qualitative agreement with the appearance in quasar absorption line spectra of a well-advanced 
state of structure formation at redshift z ~ 5, but the large amplitude of the density fluctuations may 
produce significant small-scale anisotropy in the CBR. The analysis of the CBR anisotropy in second 



order perturbation theory was introduced by Ostriker & Vishniac (1986). Vishniac's (1987) more detailed 
investigation has been confirmed by Hu et al. (1994), Dodclson & Jubas (1995), and [Hu fc White (1997) . 
The analysis and numerical evaluation are extended in Persi (1995) and Persi et al. (1995). The application 



of perturbation theory may be uncertain if early structure formation produced highly nonlinear departures 
from homogeneity at the epoch of last scattering, however. To investigate this we have developed a 
nonperturbative model for scattering in a cloudy distribution of plasma. The expression for the CBR 
anisotropy ST/T in the strongly cloudy limit bears a close resemblance to the Ostriker- Vishniac relation, 
and the numerical results for 5T/T accordingly are similar. 



Our analysis, which extends previous discussions by Hogan & Partridge (1989), Peebles (1990) 



Aghanim et al. (1996)| and pruzinov k Hu (1998) , assumes the CBR was last scattered by free electrons 



in clouds with density contrast well above unity, so the mean free distance t f between intersections of gas 
clouds along a line of sight is much larger than the typical cloud size d c i- The simplification offered by this 
limit is that the details of the matter distribution and motion on scales from d c \ to tf are not important, 
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because a line of sight on average samples only one cloud over the distance tf. Thus we can model the 
clouds as an inhomogeneous random Poisson process. In the process the joint distribution of cloud motions 
and scattering optical depths as a function of position along the line of sight is determined by the plasma 
density and velocity fields averaged through a window of width tf. If tf is larger than the scale of non- linear 
density fluctuations then the CBR anisotropy 5T/T is the sum in quadrature of a perturbative contribution 
and a shot noise term. 

The next section shows the relation between the shot noise term in ST/T and the Ostriker-Vishniac 
(1986) effect. The model of the clouds as an inhomogeneous random Poisson process is presented in §3, 
and a simplified treatment of the effect of correlated cloud motions is discussed in §4. In §5 we present 
numerical examples of the expected CBR anisotropy. 



2. Shot noise and the Ostriker-Vishniac anisotropy 

I I 

Ostriker & Vishniac (1986) showed that the CBR anisotropy at angular resolution 9 produced by the 

peculiar motion of the plasma at last scattering at proper expansion time t s is 

Here rg is the proper length subtended at the epoch of last scattering by the observed angle 9, and Sg 
and vg are the rms values of the density contrast and the line of sight component of the peculiar velocity 
averaged through a window of width rg . Ostriker & Vishniac based this relation on the lowest nontrivial 
order of perturbation theory. When the CBR is scattered in discrete and well-isolated clouds of plasma the 
analogous expression is obtained as follows. 

Suppose the CBR is scattered in clouds with typical optical depth 

r = <r£ e , (2) 

where a is the Thomson cross section and E e is the characteristic free electron column density of a cloud. 
If r ;$ 1 the CBR temperature perturbation observed at angular resolution smaller than a cloud size may 
be approximated as 

(3) 

In the integral n e and v are the free electron number density and streaming velocity as functions of position 
along a line of sight, and the integral is over a proper distance comparable to the expansion time t s at the 
epoch where the universe last is optically thick, an e t s ~ 1. (Here and below the velocity of light set to 
unity.) The sum is over N ~ r^ 1 ~ t s /tf clouds, where tf is the mean distance between cloud intersections 
along a line of sight. If tf is larger than the nonlinear clustering length then the sum in equation adds as a 
random walk, and the mean square temperature perturbation along a line of sight is 

(4£) 2 ~ w V 2 iV ~ v 2 t. (4) 

If the CBR is observed in a beam with angular size 9 larger than the angular diameter 9 C \ characteristic of 
clouds with proper width d c \ the mean square anisotropy is 
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because we are ignoring correlations among cloud positions and motions. As in equation (Q), r§ is the 
proper length subtended by the beam size 6 at last scattering. The number of free electrons in a cloud is 
~ S e d^j, so the number density n c \ of clouds satisfies 

E e d^n c i ~ n e ~ l/{t s a) , (6) 

where n e is the large-scale mean number density of free electrons. Equations (||) and @ give 

This is the same form as the Ostriker-Vishniac expression in equation (Q), where v is the rms line-of-sight 
peculiar streaming velocity of the material within a cloud and Sg is the shot noise contribution to the rms 
value of the fractional fluctuation in the number of clouds found within a randomly placed sphere of width 
re- 
in the next section we generalize this analysis to an inhomogeneous random Poisson process that takes 
account of the large-scale correlations in cloud positions and motions and the cosmic evolution of the cloud 



parameters. The former have little effect on the small-scale CBR anisotropy, for punyaev ( 1978)| noted that 
the CBR averages the perturbations of the peculiar motions across the Hubble length t s at the epoch of last 
scattering, suppressing the contribution to the anisotropy. He also showed that the contribution is further 
suppressed by the anticorrelation of the peculiar motions on either side of a density fluctuation. This effect 
was later rediscovered and confirmed by Kaiser (1984)[ 



3. Poisson model for the scattering clouds 

3.1. The inhomogeneous Poisson process 

A suitable approximation to the transfer equation for the CBR temperature perturbation observed at 
epoch t along a line of sight through a given electron distribution n e and streaming velocity field with a 
line-of-sight component v is 

4j£ = J an e vdtexp — J an e dt' . (8) 

This expression ignores the gravitational perturbation to the CBR, which is small in the small angle limit 
of interest here, and the inhomogeneity in the space distribution of the radiation at last scattering, for the 
observations tell us 5T/T ~ 10~ 5 at last scattering, while the velocity term in equation (g) is v/c ~ 10~ 3 . ft 



also ignores the relativistic correction of order ~ (v/c) 2 discussed by Hu et al. (1994) and Dodelson fc Juba£ 



(1995) , for in the small and relatively cool clouds to be expected at high redshift this correction is on the 
order of v/c ~ 10 -3 times the Ostriker-Vishniac term. (The (v/c) 2 term is induced by the velocity- velocity 
nonlinear coupling, while the perturbative analog of the above equation comes from the velocity-density 
coupling). 

We are assuming the CBR is last scattered in well-separated clouds of plasma. If the optical depth of a 
cloud on the line of sight is r, the probability the cloud has scattered a particular photon into the beam is 
1 — e~ r . If the peculiar velocity of the cloud along the line of sight is v the contribution to the integral over 
t in equation (||) by this cloud is 
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whcre the integral runs over all the plasma intersected by the line of sight subsequent to this cloud. Thus 
in our cloud model we rewrite equation (0) as 



4j£ = ^2n iab (l - e r ")v 6 exp- ^ n kcd T c , (10) 

6 k>i,cd 



where Uiab = 1 if there is a cloud with optical depth in the range r a to t q + dr a and peculiar velocity in the 
range Vf, to Vb + dv b at position U to i$ + dii along the line of sight, and rii a fc = otherwise. 

Equation ( |l0| ) represents the strongly nonlinear cloud-like distribution of the scattering plasma by the 
inhomogeneous random Poisson process rii a b. The inhomogeneity of the process includes the time evolution 
of the mean cloud properties and the large-scale structure in the distribution and motion of the clouds. 
Thus we compute statistical averages in two steps: first average over the Poisson process for given smoothed 
fields that represent the large-scale structure, and then use perturbation theory to compute the average 
over the ensemble of smoothed fields. 

For a given realization of the smoothed fields the expectation value of the Poisson process in 
equation (|l0|) is 

(riiab) = fiabdr a dv b dti/tf(i) . (11) 

The probability of finding a cloud at the position labeled by in the interval dti is dU/tf, where tf(i) 
is the local mean free distance between cloud intersections in the given large-scale structure. If there is 
a cloud at U the probability distribution in its optical depth and velocity is fi a b, with the normalization 
J fiabdr a dvb = 1. Thus the average over the Poisson process of the optical depth of a cloud found at point 
ti on the line of sight is 

= / fiabT a dT a dvb = an e (i)t f (i) . (12) 



In the last expression the mean plasma density n e (i) is a function of position t% along the line of sight, 
and averaged over the mean free path tf(i). If the mean optical depth (r), is independent of position this 
equation just says the smoothed electron density varies inversely as the mean free distance between clouds. 

The mean value of the product of the electron density and streaming velocity v(i) at U is 

(n e v) = -^-—'S^(n ia b)T a Vb 
adti 

ab 



l 



at f (i) 



dT a dv b T a V b fiab 



= nS) • (13) 

The last line, which uses equation (fl2"|), expresses the average of n e v across the Poisson process as the 
product of the smoothed electron density field and the smoothed cloud velocity field weighted by the cloud 
optical depth. 



3.2. The variance of the CBR temperature 



In the computation of the mean of (8T/T) 2 we consider first the average over the Poisson process for a 
given realization of the smoothed fields that describe the large-scale structure of the plasma distribution. 
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Because riiab m equation (JO]) is either zero or unity the mean of any positive power of riiab is the same 
as the mean of riiab- By expanding the exponential as a power series, one sees that 

(e- ni " tT °) = l-(n lab )(l-e- T "). (14) 

Averages across the Poisson process may be computed independently for each different cell label iab, so the 
expectation value of the visibility function in equation ( |Io| ) is 



rah 

= H[l-(n lab )(l-e-^)] 

= exp- [dU(l- (e- r )i)/t f (i). (15) 



In the last line, the term e~ Ta is averaged over the distribution /j a & of optical depths for clouds found at 
position i along the line of sight, as in equation (fi1|). If the cloud optical depths are large, so (e~ T ) <C 1, 
equation (15) is just the probability the line of sight intersects no clouds in the Poisson process. If (r) <C 1 
equation (15) is the usual expression in perturbation theory, 

Here n e is the local smoothed free electron number density defined in equation (|l2|). 



The average of the transfer equation ( |10|) across the Poisson distribution for a given realization of the 
large-scale structure is 



ST 
T 



'" ((l-e-» t exp-/ t0 -^-(l-(e-». (17) 



o */(*) u Jt */(*') 



Here again the averages are over the distribution f ia b of optical depths and line of sight velocities of clouds 
found at a given position t along the line of sight. If the cloud optical depths are small, we can use the 
expression for the smoothed electron number density n e (t) in equation ( |l2| ) to rewrite equation ( [i"7j ) as 

= f° on e dtv{t)e-$> n < dt ' , v(t) = (rv) t /(r} t . (18) 



o 

This is the form one would write down in perturbation theory, where the smoothed field v(t) as a function 
of position t along the line of sight is the mean peculiar velocity weighted by the cloud optical depth, as in 
equation (|l3|). 

In the expression for the mean square value of the temperature perturbation along a line of sight we 
have to consider separately the squared terms and the cross terms from the sum over the space position 
index i in equation (|To|) . Thus we write the average over the Poisson process as 

((f) 2 ) = mi+mi 

As we now discuss, the squared terms generalize the shot noise in equation (^) and the cross terms 
approximate the variance of 5T/T in perturbation theory for the smoothed fields. 



The sum of the squared terms in equation ( |Io| ) is 



(f) 2 = E ("namcdO- - ^"X 1 - z~ Tc >bv d J] e - 2n ^ T '} . (20) 

iabcd k>i,ef 
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In the last product we have taken account of the fact that each exponential factor e~ nk " fT ' ! appears twice, 
once from each factor of 5T/T. As in equation we can compute separately the average over the Poisson 
process for each different cell iab. Also, we have 

(niabnicd) — {n ia b)S a cSbd , (21) 
because the terms where a ^ c or b ^ d are of second order in dti . Thus equation (|20|) is 



<4?) - - [°T;^ 1 -°- T ^-Cm (1 - {e ~ 2T) ' ) ■ (22 > 

We complete the discussion of this shot noise term in §3.3, after dealing with the cross terms in equation jl9|). 
The sum of the cross terms in equation ( JTo| ) is 

(f)c = 2]T£(l-e-^)(l- e -^H« d (n Mb } 

i<j abed 

x (n jcd l[e- n ^) [] \{(e- n ^) JJ (e~ 2n <- T ~) ■ (23) 

ef i<k<j gh l>j,mn 

The first expectation value in this equation refers to the earlier time U along the line of sight. The second 
expectation value contains the exponential factors at the later time tj along the line of sight. In the next 
group of expectation values the exponential factors are evaluated at times between ti and tj . In the last 
group, at times greater than tj, each exponential factor appears twice, from the two factors of ST/T. We 
have 

(n jcd J] e -»W) = (n jcd e- n ^ T °) = (n Jcd )e- r ' . (24) 

ef 

The first step follows because the exponential factors with ef not equal to cd introduce terms of second 
order in dtj . The last step follows by the argument used in equation (|lj) . With equations ([ll]) and ([l5]) we 
get 

^ 2 = r dtdt 
KT)c Ji<j tf(i)t f (j) 

I = (((l-e-nvMl-e-^e-^j), 

J = e -J!<*&-i°-*)V*f e -f l>j dt « 1 -l°~' T M t f . (25) 

In this equation the outermost brackets, (...), denote an average over an ensemble of large-scale velocity 
fields, while the inner brackets, (. . denote an average over the Poisson process for each cell iab for a 
given realization of the ensemble of velocity fields. 

As in equation (|lj), when the cloud optical depths r a are much less than unity we can rewrite 
equation (E5|) as 



^ - I'' dmja'n^U) <<7>'<™>J> e -J k>i ™^ e -Li™° dt '. (26) 



Tie / "«'i u >»3 



\Vi\ T /3 



This is the expression one would write down in perturbation theory based on the smoothed velocity and 
density fields defined in equations (|l^), (|l3|), and (|l8|). If tf is large compared to the scale of nonlinear 
clustering, this smoothed velocity field is well approximated by linear perturbation theory. 
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3.3. The shot noise contribution to 8T/T 

We consider here some simplified forms for the shot noise term in equation (|2^). We assume the 
universe is ionized back to high redshift where the scattering optical depth to the present is large. If the 
distributions of r a and Vb are independent of time, the integral is 

[t) s - 1_< e -2r) • V { ) 

We see that if r ^> 1 the temperature perturbation is just the rms velocity of the last cloud along the line 
of sight. If the dispersion in r is small, equation (|27|) is 



(4£) 2 = (u 2 ) L^I = (u 2 ) ^_ (2g) 

\TJ S \ i\ + e -T \ '2-q x ' 
where the probability a photon is scattered by a cloud is 

q = l-e- T . (29) 



Another derivation of equation (p8|) is in Peebles (1990) 



In the numerical examples in §5 the cloud optical depths r are small. Here, following equations ( |l8| ) 
and (|2g), we can write equation (|22] ) as 

4£) 2 = f dt ^ - ^ an e exp -2 / an e dt' . (30) 



V T 

JO V'/ Jt 

This is the average across the Poisson process for a given realization of the smoothed fields that represent 
the large-scale structure. In linear perturbation theory the average across the smoothed fields just replaces 
the density and velocity dispersion with the global mean values as functions of world time. One could 
numerically evaluate the resulting time integrals, but if {t 2 v 2 } j (t) is not a rapidly varying function of time 
the integral is well approximated by the value of this expression at the epoch of last scattering. In this 
approximation the shot noise term is 

The factors in the last expression are suitably weighted values of the mean cloud optical depth and mean 
square peculiar velocity evaluated at the epoch of last scattering of the CBR, under the assumption that 
the baryons are concentrated in clouds back to this epoch. 

To describe the effect of scattering at lower redshifts we use the mean scattering optical depth, 

f(z) — / ah e (dt/dz)dz, (32) 
Jo 

where n e is the cosmic mean density of free electrons in optically thin clouds. In the "saturated" case, the 
integral in equation ( |32| ) increases indefinitely with z, and the redshift of last scattering, z s , is defined by 
f(z s ) = 1. Equation (pTj) assumes this saturated ionization case. Another possibility is that f(z) does not 
grow indefinitely with z but rather reaches a maximum value f < 1 at z = z*, and then remains constant 
back to the hydrogen recombination epoch. Here the approximation to the shot noise term in equation ( |30| ) 
changes from equation ( |3l| ) to 

(^)! = ^(Irf • (33) 
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The plasma streaming velocity v s at last scattering, at redshift ~ z*, has been scaled to the present value 
v by the factor v s /v a computed in linear perturbation theory. It will be recalled that the typical optical 
depth per cloud is r s . 

To compare our results with experimental data and theoretical predictions for primary anisotropies, 
it is convenient to expand the CBR temperature two-point angular correlation function in Legendre 
polynomials, 

c{@) = (f (i)f (2)) = J2 ^r c ? p t ( cos e )< ( 34 ) 

where is the angle between directions (1) and (2), while 

Ci = (Kf > , (35) 

and a™ is a coefficient in the spherical harmonic expansion of the CBR temperature distribution. The 
shot noise contribution to the CBR temperature autocorrelation function approximates a step function, 
£7(0) — (ST/T) 2 at angular separation smaller than the typical angular size 9 C \ of a cloud at redshift z s , 
and C(0) = at 9 > 9 C \. The angle 9 C \ subtended by a cloud of diameter d c \ at redshift z s is 

0ci = H d d (l + z s )/y, (36) 

where y is the usual dimensionless angular size distance (Peebles 1993, eq. [13.29]), 

y(z s ,VL,K,H ) = H a r , (37) 

with r equal to the comoving angular size distance to a cloud at z s . To avoid ringing in ^-space, we replace 
the step function-like C(8) by a Gaussian, 

£7(0) = (ST/T)l cx P (-0 2 /e c 2 1 ) . (38) 

If the shot noise dominates, and £ ^> 1, the resulting power spectrum of CBR fluctuations is given by the 
usual Hankel transform, 



POO 

Ci = (^)l / exp(-e 2 /9 2 cl ) J (0^)2^0 dQ = n8 2 cl (4£)' | exp( 
Jo 



i 2 2 cl /4) , (39) 



where J is the zeroth order Bessel function. In these models the power per octave, ~ £(£ + l)Ce, rises 
with a growing multipole number like £ 2 , reaching its peak value « (5T/T) 2 near I « 0~, , and then drops. 
This maximum signal is diluted in experiments with antenna beamwidth greater than 9 C \. The mean square 
CBR temperature anisotropy averaged through a Gaussian beam response of dispersion 9/2, corresponding 
to full width at half maximum of 1.29, is 

(ST) 2 = c\ (ST) 2 (401 

V T IB Q2 + Q2^ \ T ) s ■ \ W ) 

The image of the microwave sky produced with a radiotelescope beam 9 3> 9 C \ will be resolution limited, with 
the rms 5T/T reduced by the factor 9 c \/9. Equations ([53]) and (^) are the analog of the Ostriker-Vishniac 
expression in equation (Fy), as we discussed in connection with equation (0). 
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3.4. The contribution from correlated motions 



Here we consider the term induced by the correlations in the peculiar velocity field. We ignore the 
correlations of cloud positions in our analysis for the following reason. 

If the power spectrum of density fluctuations at zero redshift is given by P(k), where k is the 
wavenumber, then the linearized continuity equation says that the present power spectrum of the peculiar 
velocity field is P v {k) = H 2 a 2 J 2 (n, A)P(k)/k 2 , where / » Q 6 (e.g. |Peebles 1993j §13). Because of 
the factor k~ 2 the coherence length of the density field is always significantly shorter than that of the 
velocity. Moreover, density correlations tend to decrease more rapidly with increasing redshift than velocity 
correlations. Thus in an Einstein-de Sitter universe, P(k, z) — P(k)(l + z)~ 2 while P v (k, z) = P v (k)(l + z)^ 1 , 
so at high redshift 

(v z /v ) 2 = P v (k,z)/P v (k) » P(k,z)/P(k). (41) 

A similar inequality holds for a wider class of models, with < 1 and/or A ^ ( Peebles 1993| , §13). For 
all models under consideration here, the left-hand side of the above inequality is larger than the right-hand 
side by a factor ~ z s 3> 1. In other words, the spatial correlations between the clouds can be safely ignored. 

As in §3.3, to simplify our calculations further, we now assume the dispersion in the cloud optical 
depths is small and that the probability of Thomson scattering per cloud is given by equation (^9|) , where q 
and r do not depend on redshift (while in more realistic models r may be a function of the cloud and the 
redshift). Instead of the cloud coordinates along the past light cone, (U, tj), it is more convenient to use the 
mean and relative position, given by 



t = t(i,j) = (U+t^/2. 



u(i,j) = U - tj 



(42) 



as in the usual derivation of the Limber equation (see e.g., Peebles 1993, §7). After this coordinate 
transformation, relations (EEa) and (K9) become 



and 



(f): 



(fl)* = 2, 2 (1-,) 



dt -(2„-„ 2 
■ e 



*/(*) 



(2q-q 2 ) J'" dt'/t f (t') 



n(o,i). 



dt -(2a- 



e 



(2?-? 2 ) f t ° dt'/t f (t') 



10 tf(t) 

Here n is the radial component of the velocity correlation tensor, 



du 



-uq/tf (*)!!(«, t) 



(43) 



(44) 



U{u,t) = (v z /v ) 2 U(r) 

u is the proper separation, r = u/a(t) is the comoving separation, and a(t) is the expansion parameter. 
The function n(r) can be expressed in terms of the power spectrum, 



(45) 



n(r) = 



2tt 2 



P(k) 



jo(kr) - 2 



ji{kr) 



kr 



dk . 



(46) 



where the ji are spherical Bessel functions ( |Groth, Juszkicwicz, fc Ostrikcr 1989 , Peebles 1993 ). At r = 0, 
the value of the expression in square brackets is 1/3. Hence, the source term of the shot noise contribution 
(i.e., the one-dimensional velocity dispersion) is 



n(o,t) 



6tt 2 



P(k)dk. 



(47) 
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In our analysis of the relative importance of the contribution from correlated motion, this expression should 
be compared with the inner integral in equation ( [44| ) , which can be written as 

^e-^m(u,t) = (V) 2 {a f° f)2 ^ P(k)W(kt f /aq)dk, (48) 
Jo tf \ v oJ 2vr 2 q 7 

where the window function W is 

/>oo 

W(k) = k- 1 / dxe~ x ' K \j (x) - 2j 1 (x)/x] . (49) 



The variable n(k,t) = ktf/aq is the ratio the photon mean free path, tf(t)/q to the proper length a(t)/k 
corresponding to the wavenumber k. The integral ( fi"9| ) can be expressed in terms of elementary functions 
(see appendix A), 

W(k) = kT 2 - KT 3 arctan(K) » (3 + k 2 ) -1 . (50) 

The last expression is exact at k — and in the limit k — > oo. For any intermediate values of k, the 
accuracy of the approximate expression is better than 20%. The window function W[n(k, t)] acts as a low 
pass filter, damping the contribution to (ST/T) C induced by correlated motions on characteristic scales that 
are shorter than the photon mean free path. Whether or not (ST/T) C can be neglected in comparison to 
(5T/T) S thus depends on the ratio of the velocity coherence length to the photon mean free path at z s . To 
be more specific, let us model the power spectrum as 

A k 

Pm = (TTHf' ,51) 



where A, r v and (3 are constant parameters. The APM data suggest (3 — 1.2 and r v = 33 h x Mpc (Baugh fc 



Gaztahaga 1996). We trade precision for simplicity and for the purpose of our order of magnitude estimates 



set 



0=1, (52) 



To make this spectrum applicable for our linear theory expression for n(r), we need a high wavenumber 
cutoff to account for the stabilization of clustering by virialization, occurring on comoving scales smaller 
than k~l -C r v , where the density contrast exceeds unity. We therefore set P(k) = for k > k n i- Now the 
integrals in equations (^) and (E§) become trivial, and the ratio of the correlated source term to the shot 
noise source term is 

J k -'PWdk |ln[(^A 2 + 3)/3(^ + l)]| 



/ fc "' Pdk q\{X/r v f-'i\ ln(l + fc2 ;r 2) 



(53) 



where A = tf/qa, and r v /X is the ratio of the velocity coherence length in proper coordinates, r v a(z), to 
the mean free path of the CBR photons, tf(z)/q. The correlated term becomes negligibly small when the 
condition 

- f^) 2 « 1 (54) 
q \tf(Zs)z s J 

is satisfied. Note that to make the correlated term dominant, it is not enough to require that 
r v a(z s ) 3> tf(z s )/q; the clouds have to be optically thin as well (q <C 1). The coherence of the velocity field 
can amplify the net CBR temperature fluctuation by causing the effects of individual clouds to add with the 
same sign and not as a random walk. This, however, can happen only if (1) the velocity coherence length 
is larger than the photon mean free path and (2) the rescattering cutoff in the sum appears sufficiently far 
away from the observer so that amplification by coherent motions is not destroyed by rescattering. The 
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second condition means good visibility, or q -C 1. Indeed, when r v /X 3> 1, the ratio of the correlated term 
to the shot noise term is 

Jo"' PWdk ^ MV3r v /\) 
J^' Pdk ~ 3q\n(k nl r v ) ' 

Hence, unless the optical depth per cloud is small, the correlated term will be subdominant, even in the 
large- r v limit. 



4. A simplified model 

We present here a simplified model that helps clarify the physics behind the contributions to ST/T 
from shot noise and correlated motions of the gas clouds. The model assumes the universe is static, each 
cloud has the same scattering probability q = 1 — e~ T (eq. |2^]), and the mean free path between clouds is 
large enough that we can neglect the correlation in cloud positions. The model takes account of the broader 
correlation of peculiar velocities. We characterize the strength and range of the velocity correlations by an 
'effective' correlation function, 

0(r) - al er r ^ , (56) 

where r v is the velocity coherence length and erj; = 11(0) is the one-dimensional velocity dispersion. The 
assumption that the correlation of cloud positions may be neglected means the probability distribution of 
relative separations r between pairs of clouds, intersected along a line of sight, and numbered i and i + j, is 
the Poisson expression 

dp(r,j) 

where 

r f = 

Jo 

is the mean separation between clouds along the line of sight, n is the spatial density of clouds and d c i is 
the cloud diameter and the cloud numbers increase with distance from the observer. The perturbation to 
the CBR along a given line of sight is (Peebles 1990) 

4£ = qv x + q(l - q)v 2 + q(l - q) 2 V 3 + ... (59) 

The line of sight component of the j th velocity is Vj, counting back from the present. The prefactor 
q(l — q)- 1 is the probability q that a photon is scattered into our line of sight by the cloud j, multiplied 
by the probability that it is not scattered out of our line of sight by the remaining j — 1 clouds in 
the foreground. If the universe is ionized back to large redshift equation ( |59| ) in effect is a convergent 
infinite series, meaning the influence of consecutive terms decreases with increasing cloud number j. The 
attenuation by rescattering out of the line of sight means the CBR perturbation is determined by the last 
~ ./V cloud intersections, where 

N ^ q- 1 . (60) 

The attenuation by rescattering is most severe when the clouds are opaque (q = 1), and all j > 1 terms 
vanish. The nearest cloud obscures the more distant ones, the anisotropy is determined only by the velocity 
of the nearest cloud, ST/T = v\, and 

((f) 2 ) (61) 



j'-i 



( r / r f) 
(j-1)! 



-r/r f 



dr 
'7 



rdp(r, 1) = , 



(57) 
(58) 



-12- 



Let us now return to the more interesting case when q < 1. Squaring the expression (p9]), and averaging, 
we get 



oo oo 

'f) 2 > = ^EEf 1 -^^),, (62) 



where (n)„ = ex 2 and for j > i the mean of the velocity autocorrelation function (eq. K^j) for the Poisson 
distribution in the distance between cloud intersections i and j (eq. p7| ) is 

mV -a 2 [°° — (—Y * 1 - M (l/r /+ l/r„) _ /g 3l 



(64) 



With this expression the sums in equation (B2| are readily evaluated to get 



((f) 2 ) = ^ 



9 



i + 2(1 - g)(r v /r f ) 
l + q(r v /r f ) 



The first term in the square brackets comes from the sum of the squared terms in equation ( p2[) and the 
second from the sum of the cross terms. 

If the velocity field is uncorrelated, r v = 0, the contribution from the cross terms vanishes and 
equation (|64| ) reduces to equation (p8|) . The cross terms also vanish if the nearest cloud along the line 
of sight is opaque, q = 1, and the anisotropy reduces to the "single cloud" limit (^lj), independent of the 
velocity correlation length. The "single cloud" limit also applies when the velocity correlation length is 
large, qr v r/, because the photon mean free path samples a single coherently moving set of clouds. And 
more generally, a positive velocity correlation reduces the effective number of statistically independent steps 
in the random walk and increases the rms temperature perturbation. 

We note finally that one can take account of the Hubble expansion by replacing equation ( p7| ) with an 
inhomogeneous Poisson process, 

di ( f° df/Ut'))*' 1 fto ,„ 
a m at \Jt i i\ jj _ / dt' t f (t') 

from which it is an interesting exercise to rederive equations (|43|) and (|44|). 



5. Numerical examples 

Since we know very little about the history of structure formation at redshifts z <L 5 a good strategy is 
to consider examples of what might have happened and how it would have affected observables such as the 
CBR temperature fluctuations. To keep our discussion of possibilities simple and definite we choose one set 
of cosmological parameters, 



n = o.2. 



ill 



0.05, 



h = 0.7 = Ho/100 km s _1 Mpc~ 



(66) 



The low value of is in line with the observational evidence ( Freedman et al. 1994 ; Peebles 1997b| 



Pcrlmutter et al. 1998; and references therein), and cosmogonies with the early structure formation we have 



assumed in this analysis seem to be most promising in a low density cosmological model (Peebles 1997a). 
To emphasize the scattering effect we have adopted a baryon density parameter near the high end of the 



- 13 - 



range now under discussion (Copi, Schramm, & Turner 1995). Following Navarro, Frenk, & White (1997) 



we assume that at any epoch matter is concentrated in clouds with density contrast 

5 ~ 200. (67) 

Then a cloud mass M fixes a characteristic cloud diameter, 

dd = [6M/(^M)] 1/3 , (68) 

where p s is the mean mass density at the scattering epoch z s . A measure of the physical length scale of 
nonlinear mass fluctuations is 

15Mpc /£ s \ 2 /( 3 +") 



h(l + z s ) \D 0/ ' ' l!>,,! 

A sphere with this diameter contains roughly the cloud mass M. We have normalized to the rms fluctuation 
of galaxy counts, 8N/N ~ 1 in a sphere of diameter 15ft. -1 Mpc. The assumption that galaxies trace mass 
agrees with our low value for the density parameter fi. Equation (|69| ) is scaled by the power law power 
spectrum 

P(xk m , (70) 

to the epoch z s using the linear perturbation theory growth factor for the density contrast, Sp/p oc D (t). 
Finally, we take the present rms peculiar velocity to be 

u o = 600 kms" 1 , (71) 

and we scale this velocity with time by the linear perturbation theory relation v oc adD/dt. 

We consider first the effect of scattering by clouds at epoch 1 + z s = 10. In a cosmologically flat model 
the density and velocity growth factors are 

l + z s = 10, D./Do = 0.14, v s /v o = 0A9. (72) 

For a mass characteristic of the central part of a large galaxy the power spectrum index (eqs. |6{J, |Q), the 
model parameters are 

Mci = 1 x lO n M 0) m = -lA, 77 =4 Mpc, 

d c i = 33 kpc, y = 2.5, 6 C \ = 6.4 arc sec . (73) 

Here rj is the mean free path for intersection of clouds, d c \ and 6 C \ are the cloud diameter and the angle 
it subtends, finally, y is the usual angle size distance parameter. (For comparison, note that for A = 
and z s ^> Or 1 , the angular size distance is y = 2/0.) The shot noise contribution to the CBR anisotropy 
(eqs. @ and @) is 

8T/T{6 = 1 arc min) = 2 x 10~ 6 . (74) 

We are assuming all baryons are in optically thin plasma, so the mean optical depth for scattering back 
to this redshift is f = 0.1; the probability of scatternig per cloud is r s = 0.004. The contribution to the 
Compton- Thompson parameter in the CBR spectrum by the motions of the galaxies is y c ~ fv 2 s = 1 x 10~ 7 , 
well within the COBE bound (Fixsen et al. 1996; note that v s here is expressed in units of the velocity of 
light; in our units c = 1). 

With the parameters in the above example the mean baryon density within a cloud is n ~ 0.06 protons 
cm~ 3 , and the net mass density is equivalent to about 0.2 protons cm~ 3 . These numbers may not be 
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unreasonable for a protogalaxy. The power law index in e quation (|7q) is comparable to that of the most 
recent version of the isocurvature CDM model, m ~ —1.8 ( Peebles 1997a ). If the characteristic cloud mass 
is reduced to M c \ = 1 x 10 10 M Q with all other parameters the same we get m = —1.8 and 8T/T = 6x 10~ 7 
in a window of one arc minute diameter. In an open universe (with zero cosmological constant) with 
M c i = 1 x lO n M0 and all other parameters the same ST/T is 1.5 times the value in equation d73). 



The relative size of the contribution to 5T/T by correlated motions is determined by equation 
With the numbers in equation ([73]) the scattering probability in a cloud is q = 0.004, so if the parameter 
q(r v / z s rf) 2 were of order unity the comoving velocity coherence length would have to be r v ~ 700 Mpc, 
much larger than that suggested by the APM power spectrum (eq. plj; Baugh & Gaztahaga 1996) or even 
the most radical interpretation of bulk flow observations. That is, 5T/T in our model is dominated by the 
shot noise. 

At a larger redshift in a cosmologically flat model, 

1 + z s = 30, D s /D = 0.047, v s /v a = 0.28, (75) 

the probability of scattering is t s = 0.007 per cloud and the mean optical depth is f = 0.6 if all baryons are 
in optically thin plasma, and a cloud model is 

M cl = 1 x 1O 9 M , m = -1.4, 77 = 300 kpc, 

d c i = 2.4 kpc, y = 3.1, 9 C \ = 1-1 arc sec, (76) 

and 

5T/T(6 = 1 arc min) = 5 x 10~ 7 . (77) 
At still larger redshift in a cosmologically flat model, 

1 + z s = 50, D S /D D = 0.028, v s /v a = 0.22, (78) 

we get t s = 0.004, f = 1, and a cloud model is 

M c i = 1 x 10 7 M Q , m = -1.7, rf = 40 kpc, 

d c i = 0.3 kpc, y = 3.3, 6 C \ = 0.2 arc sec , (79) 

and the contribution to the CBR anisotropy is still smaller, 

ST/T(9 = 1 arc min) = 7x 10~ 8 . (80) 

Here the contribution to the Compton-Thompson parameter is y c ~ 1 x 10~ 7 . 



The temperature fluctuations in these examples are well below the measured bounds (Readhcad ct al 



198S ; Fomalont et al. 1993| ; [Bubrahmanian et al. 1995 ; Church ct al. 1997 ; Partridge et al. 1997 ; Andreani 



1994 ). In Figure 1 we plot power spectra for CBR anisotropies for the above three models, calculated from 



equation (|39|) . One sees that in the models we consider measurements that resolve the clouds should detect 
secondary temperature fluctuations as large as primary anisotropies at their maximum. 



6. Discussion 



The simplifying assumption for this analysis is that structure formation is so well advanced at the 
epoch of last scattering of the CBR by free electrons that the mean distance t f between intersections of 
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clouds along a line of sight is large compared to the scale of nonlinear mass fluctuations. This allows us 
to model the clouds as an inhomogeneous random Poisson process determined by the mass density and 
peculiar velocity fields smoothed through a window of width 1 1 , and it leads to the shot noise contribution 
to the small-scale CBR anisotropy in equation (p3|). This approach is motivated by the isocurvature CDM 



model for structure formation (Peebles 1997a), in which structure formation could commence at decoupling 
at redshift z ~ 1000. A second important motivation has been to complement the usual perturbative 
analysis of the effect of the nonlinear growth of small-scale structure. The similarity of results from the 



perturbative ( Persi et al. 1995 ) and nonperturbative approaches leads us to believe we have reliable 



methods for estimating the effect of early nonlinear structure formation on the CBR anisotropy. 

The observations of young galaxies and the intergalactic medium at z ~ 3 indicate a situation 
intermediate between the perturbative and nonperturbative cases. The damped Lyman-a systems contain 
a significant baryon fraction, and the mean distance between intersections of these clouds is large (at z = 3 
it is longer than the Hubble length). There also is a significant baryon fraction in the Lyman-a forest, 
and these clouds have a relatively short mean free distance. It is not unreasonable to speculate that the 
situation at much larger redshifts similarly calls for a combination of the two approaches to the analysis of 
the angular distribution of the CBR. 

The CBR anisotropy produced by the Sunyaev-Zel'dovich (1970) effect of the hot electrons in clusters 
of galaxies, which certainly is dominated by the shot noise term, offers an important constraint on the 
epoch of collection of the intracluster plasma. The evidence from the analysis of Persi et al. (1995) is 
that this constraint does not yet rule out the early structure formation picture. And our conclusion from 
the numerical examples in §5 is that within presently known observational constraints structure formation 
could have commenced when the universe was optically thick to scattering of the CBR. 
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A. The Thomson scattering window function 



The integral ([49] ) can be expressed in terms of two hypergeometric functions ( Ryzhik fc Gradshtcyn 



1994 hereafter GR, eq. [6.621.1]), 

W(k) = {l-K 2 )- 1 ' 2 



F ( I 1.3. k 2 \ 2 p ( 1 

r \2 ' 2 ' 2 ' 1+K 2 ) 3 \2 



3. 5. k 

2 ' 2 ' 2 ' 1+k' 2 



(Al) 



One of the hypergeometric functions turns out to be an elementary function in disguise (GR, eq. [9.121.26]): 



F 



■ 1 1 . 3 . ^2 



z 2 ) = arcsin(z)/z , 



(A2) 



while the other can also be expressed in terms of elementary functions by differentiating equation ( |A2| ), and 

(A3) 



using the standard recurrence relation (Abramowitz fc Stegun 1972, eq. [15.2.7]) 



I- [(1 - z) a F (a, b;c;z)} = (1 _ z )^ F ( a + 1, b; c+ 1; z) 

oz c 
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for (a, b, c) = (1/2, 1/2, 3/2). The final result of all this is the expression (j50[). 
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Figure 1 

CBR angular power spectra, calculated from equations ([33|) and ( |39| ) for three sets of model parameters, 
described in §5. The anisotropy is induced by scattering in moving clouds of ionized hydrogen. To identify 
each of the three models, we label the power spectra with the appropriate values of the cloud angular sizes, 



